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VARIATIONAL RELEVANCE VECTOR MACHINE 



FIELD OF THE INVENTION 

This invention relates generally to data modeling and analysis, and more 
particularly to a variational relevance vector machine for such data modeling and 
analysis. 

BACKGROUND OF THE INVENTION 

Data modeling has become an important tool in solving complex and large real- 
world computerizable problems. Applications of data modeling include data 
compression, density estimation and data visualization. A data modeling technique used 
for these and other applications is probabilistic modeling. It has proven to be a popular 
technique for data modeling applications such as speech recognition, vision, handwriting 
recognition, information retrieval and intelligent interfaces. One framework for 
developing such applications involves the representation of probability distributions as 
directed acyclic graphs, which are also known as Bayesian networks, belief networks, 
and probabilistic independence networks, among other terms. 

In probabilistic modeling, usually a training data set is given that includes input 

vectors {x„}^ =1 along with a set of corresponding targets {^j^ , the latter of which can 

be real values, in the case of regression analysis, or class labels, in the case of 
classification analysis. From this training set, a model of p(t\x) is attempted to be 
inferred, with the object of making accurate predictions of t for new, unlabelled, 
examples of x. Generally, the principal challenge is to find the appropriate complexity of 
this model. Scoring alternative models by training set accuracy alone is usually 
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undesirable, since increasing the model complexity, while reducing the training set error, 
can easily lead to over-fitting and poor generalization. A more robust approach is to 
introduce a prior distribution over models, which is used in conjunction with the 
information supplied by the training data to infer the prediction model. This prior 
5 distribution, also referred to as a prior, can be explicit, such as in a Bayesian framework, 
or can be implicit in other approaches. 

One method for classification, that has also been extended to regression, is known 
as the support vector machine (SVM). Although it does not estimate p(t\x) 9 it makes 
predictions based on a discriminant function of the form 

N 

10 y (x) = ^w n K(x,x H ) + w 09 

where {w n J are the model weights and K(%>) is a kernel function. A feature of the SVM is 
that its cost function attempts to minimize the number of errors made on the training set 
while simultaneously maximizing the margin between the two classes, in the feature 
space implicitly defined by the kernel. This maximum-margin principle is an appealing 

15 prior for classification, and ultimately drives many of the weights to zero, resulting in a 
sparse kernel classifier where the non-zero weights are associated with x n that are either 
on the margin or lie on the wrong side of it. Model complexity is thus constrained such 
that only these support vectors determine the decision function. In practice, in addition to 
fitting the model to the training data, it is also necessary to estimate the parameter 

20 (usually, denoted Q which regulate the trade-off between the training errors and size of 
margin, which may entail additional cross-validation. 

A disadvantage with the SVM as a general matter is that it utilizes many kernel 
functions, and may not yield as optimal test performance as may be desired. 



Furthermore, the SVM utilizes parameters (i.e., those denoted Q, which add unwanted 
complications to the model. To address these concerns, the copending and coassigned 
patent application entitled "Relevance Vector Machine," filed on September 4, 1999, and 
assigned serial no. 09/391,093, describes a Relevance Vector Machine (RVM) that 
5 utilizes a functional form that is equivalent to the SVM, but which is a probabilistic 
model. It achieves comparable recognition accuracy to the SVM, but advantageously 
provides a full predictive distribution, and requires substantially fewer kernel functions. 
As described in this prior application, the RVM relied on the use of type II maximum 
likelihood, referred to as the evidence framework, to generate point estimates of the 
10 hyperparameters that govern model sparsity. However, because analysts desire to have 
different approaches, techniques and tools to solve a given model, there is a motivation 
for the present invention. Furthermore, the approach described here provides a closer 
approximation to a fully Bayesian treatment than has been possible previously, and this is 
expected to be advantageous for problems involving data sets of limited size. 

15 

SUMMARY OF THE INVENTION 

The invention relates to a variational relevance vector machine (RVM). The 
RVM is a probabilistic basis model of the same functional form of the SVM. Sparsity is 
achieved through a Bayesian treatment, where a prior is introduced over the weights 
20 governed by a set of what are referred to as hyperparameters - one such hyperparameter 
associated with each weight, whose most probable values are iteratively estimated from 
the data. The posterior distribution of many of the weights is sharply peaked around 
zero, in practice. The variational RVM utilizes a variational approach to solve the model, 
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in particular using product approximations to obtain the posterior distribution. In 
particular, the product approximation used is the distribution of the hyperparameters, 
times the distribution of the weights, times the distribution of predetermined additional 
parameters, that differ as to whether a classification-oriented (discrete) model is being 
5 provided, or a regression-oriented (continuous) model is being provided. For the latter, 
the predetermined additional parameters account for noise, while for the former, the 
predetermined additional parameters account for a lower bound. 

In one embodiment, a computer-implemented method includes selecting an initial 
set of hyperparameters, and iteratively updating the distribution of the set of weights, the 

10 distribution of the set of hyperparameters, and the distribution of the set of predetermined 
additional parameters, until a predetermined convergence criterion has been reach. Thus, 
the product of these distributions, as they have been iteratively updated until the 
predetermined convergence criterion has been reached, approximates the posterior 
distribution for modeling of a data set. 

15 RVM has advantages not found in prior art approaches such as SVM, As 

compared to SVM, for example, the non-zero weights in the RVM have been seen to not 
be associated with examples close to the decision boundary, but rather appear to represent 
more prototypical examples of classes. These examples are termed relevance vectors. 
Generally, the trained RVM utilizes many fewer basis functions than the corresponding 

20 SVM, and often exhibits superior test performance. Furthermore, no additional 
validation of parameters (such as C) is necessary to specify the model, save those 
associated with the basis. 
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The invention includes computer-implemented methods, machine-readable media, 
computerized systems, and computers of varying scopes. Other aspects, embodiments 
and advantages of the invention, beyond those described here, will become apparent by 
reading the detailed description and with reference to the drawings. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a diagram of an operating environment in conjunction with which 
embodiments of the invention can be practiced; 

FIG. 2 is a diagram of a predictional model; 

FIG. 3 is a diagram of a learning machine; and, 
10 FIG. 4 is a flowchart of methods according to varying embodiments of the 

invention. 

DETAILED DESCRIPTION OF THE INVENTION 

In the following detailed description of exemplary embodiments of the invention, 
reference is made to the accompanying drawings which form a part hereof, and in which 

15 is shown by way of illustration specific exemplary embodiments in which the invention 
may be practiced. These embodiments are described in sufficient detail to enable those 
skilled in the art to practice the invention, and it is to be understood that other 
embodiments may be utilized and that logical, mechanical, electrical and other changes 
may be made without departing from the spirit or scope of the present invention. The 

20 following detailed description is, therefore, not to be taken in a limiting sense, and the 
scope of the present invention is defined only by the appended claims. 
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Some portions of the detailed descriptions which follow are presented in terms of 
algorithms and symbolic representations of operations on data bits within a computer 
memory. These algorithmic descriptions and representations are the means used by those 
skilled in the data processing arts to most effectively convey the substance of their work 
5 to others skilled in the art. An algorithm is here, and generally, conceived to be a self- 
consistent sequence of steps leading to a desired result. The steps are those requiring 
physical manipulations of physical quantities. Usually, though not necessarily, these 
quantities take the form of electrical or magnetic signals capable of being stored, 
transferred, combined, compared, and otherwise manipulated. 

10 It has proven convenient at times, principally for reasons of common usage, to 

refer to these signals as bits, values, elements, symbols, characters, terms, numbers, or the 
like. It should be borne in mind, however, that all of these and similar terms are to be 
associated with the appropriate physical quantities and are merely convenient labels 
applied to these quantities. Unless specifically stated otherwise as apparent from the 

15 following discussions, it is appreciated that throughout the present invention, discussions 
utilizing terms such as processing or computing or calculating or determining or 
displaying or the like, refer to the action and processes of a computer system, or similar 
electronic computing device, that manipulates and transforms data represented as 
physical (electronic) quantities within the computer system's registers and memories into 

20 other data similarly represented as physical quantities within the computer system 

memories or registers or other such information storage, transmission or display devices. 
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Operating Environment 

Referring to FIG. 1, a diagram of the hardware and operating environment in 
conjunction with which embodiments of the invention may be practiced is shown. The 
description of FIG. 1 is intended to provide a brief, general description of suitable 

5 computer hardware and a suitable computing environment in conjunction with which the 
invention may be implemented. Although not required, the invention is described in the 
general context of computer-executable instructions, such as program modules, being 
executed by a computer, such as a personal computer. Generally, program modules 
include routines, programs, objects, components, data structures, etc., that perform 

1 0 particular tasks or implement particular abstract data types. 

Moreover, those skilled in the art will appreciate that the invention may be 
practiced with other computer system configurations, including hand-held devices, 
multiprocessor systems, microprocessor-based or programmable consumer electronics, 
network PC's, minicomputers, mainframe computers, and the like. The invention may 

1 5 also be practiced in distributed computing environments where tasks are performed by 
remote processing devices that are linked through a communications network. In a 
distributed computing environment, program modules may be located in both local and 
remote memory storage devices. 

The exemplary hardware and operating environment of FIG. 1 for implementing 

20 the invention includes a general purpose computing device in the form of a computer 20, 
including a processing unit 21, a system memory 22, and a system bus 23 that operatively 
couples various system components include the system memory to the processing unit 21. 
There may be only one or there may be more than one processing unit 21, such that the 
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processor of computer 20 comprises a single central-processing unit (CPU), or a plurality 
of processing units, commonly referred to as a parallel processing environment. The 
computer 20 may be a conventional computer, a distributed computer, or any other type 
of computer; the invention is not so limited. 

The system bus 23 may be any of several types of bus structures including a 
memory bus or memory controller, a peripheral bus, and a local bus using any of a 
variety of bus architectures. The system memory may also be referred to as simply the 
memory, and includes read only memory (ROM) 24 and random access memory (RAM) 
25. A basic input/output system (BIOS) 26, containing the basic routines that help to 
transfer information between elements within the computer 20, such as during start-up, is 
stored in ROM 24. The computer 20 further includes a hard disk drive 27 for reading 
from and writing to a hard disk, not shown, a magnetic disk drive 28 for reading from or 
writing to a removable magnetic disk 29, and an optical disk drive 30 for reading from or 
writing to a removable optical disk 31 such as a CD ROM or other optical media. 

The hard disk drive 27, magnetic disk drive 28, and optical disk drive 30 are 
connected to the system bus 23 by a hard disk drive interface 32, a magnetic disk drive 
interface 33, and an optical disk drive interface 34, respectively. The drives and their 
associated computer-readable media provide nonvolatile storage of computer-readable 
instructions, data structures, program modules and other data for the computer 20. It 
should be appreciated by those skilled in the art that any type of computer-readable media 
which can store data that is accessible by a computer, such as magnetic cassettes, flash 
memory cards, digital video disks, Bernoulli cartridges, random access memories 



(RAMs), read only memories (ROMs), and the like, may be used in the exemplary 
operating environment. 

A number of program modules may be stored on the hard disk, magnetic disk 29, 
optical disk 31, ROM 24, or RAM 25, including an operating system 35, one or more 

5 application programs 36, other program modules 37, and program data 38. A user may 
enter commands and information into the personal computer 20 through input devices 
such as a keyboard 40 and pointing device 42. Other input devices (not shown) may 
include a microphone, joystick, game pad, satellite dish, scanner, or the like. These and 
other input devices are often connected to the processing unit 21 through a serial port 

1 0 interface 46 that is coupled to the system bus, but may be connected by other interfaces, 
such as a parallel port, game port, or a universal serial bus (USB). A monitor 47 or other 
type of display device is also connected to the system bus 23 via an interface, such as a 
video adapter 48. In addition to the monitor, computers typically include other peripheral 
output devices (not shown), such as speakers and printers. 

1 5 The computer 20 may operate in a networked environment using logical 

connections to one or more remote computers, such as remote computer 49. These 
logical connections are achieved by a communication device coupled to or a part of the 
computer 20; the invention is not limited to a particular type of communications device. 
The remote computer 49 may be another computer, a server, a router, a network PC, a 

20 client, a peer device or other common network node, and typically includes many or all 
of the elements described above relative to the computer 20, although only a memory 
storage device 50 has been illustrated in FIG. 1. The logical connections depicted in FIG. 
1 include a local-area network (LAN) 51 and a wide-area network (WAN) 52. Such 
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networking environments are commonplace in office networks, enterprise- wide computer 
networks, intranets and the Internet, which are all types of networks. 

When used in a LAN-networking environment, the computer 20 is connected to 
the local network 51 through a network interface or adapter 53, which is one type of 
5 communications device. When used in a WAN-networking environment, the computer 
20 typically includes a modem 54, a type of communications device, or any other type of 
communications device for establishing communications over the wide area network 52, 
such as the Internet. The modem 54, which may be internal or external, is connected to 
the system bus 23 via the serial port interface 46. In a networked environment, program 
10 modules depicted relative to the personal computer 20, or portions thereof, may be stored 
in the remote memory storage device. It is appreciated that the network connections 
shown are exemplary and other means of and communications devices for establishing a 
communications link between the computers may be used. 

Generalized Linear Prediction Models 

15 In this section of the detailed description, generalized linear prediction models are 

overviewed, as an introduction to the relevance vector machines described in the 
succeeding sections of the detailed description. Referring first to FIG. 2, a diagram of a 
prediction model is shown. The model 200 can be a linear prediction model, which 
accepts an input 202 and provides an output 204 based on the input 202. Thus, the model 

20 200 must be constructed, which in one embodiment entails generating a prior distribution 
for the model based on a training set of input and output data input into a learning 
machine, such as a relevance vector machine (RVM) according to embodiments of the 
invention as described in succeeding sections of the detailed description. 
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Referring next to FIG. 3, a diagram for a learning machine is shown. The 
learning machine 300 accepts a training set of data x 302, and ultimately outputs a 
posterior distribution p(t\x) 304, as has been described in the background section. The 
posterior distribution 304 is determined more specifically based on an output y(x) of the 
5 learning machine 300. In a generalized linear prediction model, the output y(x) of a 
learning machine (i.e., the machine used to obtain a model that can be utilized for 
predictional purposes) can be written as a link function g{-} applied to a linearly 
weighted sum, which can be denoted as s(x), of some, usually fixed, basis functions $x), 
where x is the vector of input variables to the machine. In other words, 
10 y(x) = g{s(x)}, 

where the weighted sum is 

M 

j(x) = J]w l 4(x), 

1=1 

and Wi are the weights associated with the M basis functions. 

Relevance Vector Machines (RVM's) Generally 

15 In this section of the detailed description, an overview of RVM's is presented. 

This overview is consistent with the treatment of RVM's provided in the copending and 
coassigned patent application entitled "Relevance Vector Machines," filed on September 
4, 1999, and assigned serial no. 09/391,093. In subsequent sections of the detailed 
description, variational RVM's according to varying embodiments of the invention are 

20 described. 

The RVM makes probabilistic predictions but still retains the predictive 
performance of the support vector machine (SVM). It preserves the sparseness property 
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of the SVM. For a wide variety of problems, it has been found to actually lead to models 
that are dramatically sparser than the corresponding SVM, while sacrificing little if any 
of the accuracy of the prediction. 

For regression problems, the RVM models the conditional distribution of the 
5 target variable, given an input vector x, as a Gaussian distribution of the form 

P(t |x,w, r) = N (t\y (x,w) , r' x ) 

where the notation N(z | m 5 S) is used to denote a multi-variate Gaussian distribution 
over z with mean m and covariance S, In the above equation, t is the inverse noise 
10 parameter, and the conditional mean^x, w) is given by 

N 

j(x,w) = £w m 4(x) = w T (|>. 

m=0 

Assuming an independent, identically distribution data set X = {x n } 9 T = {tj , the 
likelihood function can be written 

15 />(r|X,w,r) = np(^|x K w ? r) 

The parameters are given a Gaussian prior of the form 

P(yv\a) = flN(w m \0,a m l ) 

m=0 

where a = {a m } is a vector of hyperparameters, with one hyperparameter a m assigned to 

each model parameter w m . In the previously filed patent application, values for these 
20 hyperparameters are estimated using the framework of type-II maximum likelihood in 
which the marginal likelihood P(T \ X,a, r) is maximized with respect to a and t. 

Evaluation of this marginal likelihood requires integration over the model parameters 
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P(T\X,a 9 r)= \p(T\X 9 ^ 9 r)P{^\a)d^ 

Since this involves the convolution of two exponential-quadratic functions, the 
integration can be performed analytically, giving 

P(T\X 9 a 9 r) = N(t\0 9 S) 

5 where t = (t\ 9 . . . , fa), and 

S = r" 1 I + 0>A" 1 O T 

in which I is the NxN unit matrix, A = diagOJ , and <t> is the Nx(N + l) design 
matrix with columns (j) m 9 so that (0) BIB = ^(x B ;x w ) . Maximization of 

P(T\X 9 a 9 r) = N(t\Q 9 S) 

10 with respect to the {a m } can be performed efficiently using an iterative re-estimation 

process obtained by setting the derivatives of the marginal log likelihood to zero. During 
the process of this optimization, many of the a m are driven to large values, so that the 
corresponding model parameters w m are effectively pruned out. The corresponding terms 
can be omitted from the trained model, with the training data vectors x n associated with 

1 5 the remaining kernel functions being termed relevance vectors. A similar re-estimation 
process is used to optimize t simultaneously with the a m parameters. 

In the classification version of the RVM, the conditional distribution of targets is 
given by 

P(t\x 9 w) = a(y) t [l-a(y)f t 
20 where <r(y) = (l + oxp(-y))~ l andy(x, w) is given by 

N 

j(x,w) = £w^ m (x) = w T (|> 
13 



Attention is confined to the case t e {0,1}. Assuming independent, identically distributed 
data, the likelihood function is obtained in the form 



in order to evaluate the marginal likelihood can no longer be performed analytically. 
Therefore, in the previous patent application, a local Gaussian approximation is used to 

10 the posterior distribution of the weights. Optimization of the hyperparameters can then 
be performed using a re-estimation framework, and alternates with re-evaluation of the 
mode of the posterior, until convergence. Thus, the previous patent application estimates 
point values for the hyperparameters. As will be described in the next sections of the 
detailed description, embodiments of the invention utilize variational inference to achieve 

15 a more complete Bayesian treatment. 

Variational Inference 

In this section of the detailed description, an overview of variational inference as 
can be performed in accordance with embodiments of the invention is given. In a general 
probabilistic model, the stochastic variables can be partitioned into those corresponding 
20 to the observed data, denoted D, and the remaining unobserved variables, denoted 0. The 



p(r|x,w)=f[^„) f "[i-^„)] 



As before, the prior over the weights takes the form 



5 




m=0 



However, the integration required by 



P(T\X,a,r)= \p(T\X,w,r)P(w\a)dv/ 
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marginal probability of the observed data, which is the model evidence, is obtained by 
integrating over 9, 

P(Z>) = \p{D,e)d& 

However, this integration typically is analytically intractable. Variational methods 
5 address this difficulty by introducing a distribution Q{0) , which, for arbitrary choice of 
Q, allows the marginal log likelihood to be decomposed into terms 

lnP{D) = L(Q) + KL(Q\\P) 

where 

10 and KL(Q || P) is the Kullback-Leibler divergence between Q{ff) and the posterior 
distribution P{0 \ D) , and is given by 

Since KL (Q \\ P) > 0 , then L(Q) is a rigorous lower bound on In P(D). Furthermore, 
since In P(D) is independent of Q, maximizing L(Q) is equivalent to minimizing 
15 KL(Q || P) , and therefore Q{0) represents an approximation to the posterior distribution 

Therefore, for a suitable choice for the Q distribution, the quantity L(Q) may be 
tractable to determine, even though the original model evidence function is not. Thus, a 
variational approach seeks to choose a suitable form for £)(#)that is sufficiently simple 
20 that the lower bound L(Q) can be evaluated and yet which is sufficiently flexible that the 
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bound is reasonably tight. Generally, a family of Q distributions is selected, and then the 
best approximation within this family is sought by maximizing the lower bound with 
respect to Q. One approach is to assume some specific parameterized functional form for 
Q and then to optimize L with respect to the parameters of the distribution. 
5 Another approach, as is described more fully in the remainder of the detailed 

description, is to consider a factorized form over the component variables in 0 9 so 

that 

The lower bound is then maximized over all possible factorial distributions by 
10 performing a free- form maximization over the Q u leading to 

exp(lnP(A<?)L 
^ A l) Jexp(lnP(A^)L^ 

where (•) denotes an expectation with respect to the distributions Q k (0 k ) for all k*i. 
If the probabilistic model is expressed as a directed acyclic graph with a node for each of 
the factors Q^) , then the solution for Qffi) depends only on the Q distributions for 

1 5 variables that are in the Markov blanket of the node i in the graph. 

It is noted that the last equation of the previous paragraph represents an implicit 
solution for the factors Q t {0 t ) , because the right-hand side depends on moments with 

respect to the Q k ^ . For conjugate conditional distributions (e.g., linear-Gaussian models 
with Gamma priors, in the case of continuous variables), this leads to standard 
20 distributions for which the required moments can be evaluated. A solution is then 
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determined iteratively by initializing the moments, and cycling through the variables, 
updating each distribution in turn using 



Controlling Model Complexity 

In this section of the detailed description, the manner by which varying 
embodiments of the invention can provide for controlling model complexity is described. 
The Relevance Vector framework generally provides a manner for solving regression and 
classification problems in which models are sought that are highly sparse, by selecting a 
subset from a larger pool of candidate kernel functions, one for each example in the 
training set. Continuous hyperparameters are used to govern model complexity, thereby 
avoiding the intractable problem of searching over an exponentially large discrete space 
of model structures. 

One manner to remove superfluous parameters it to use a pruning prior given by a 
Laplacian of the form 



However, such a choice of prior does not lead to a tractable variational treatment, since 
the corresponding variational solution cannot be evaluated analytically. 

Therefore, an alternative framework is used, based on a hierarchical prior of the 

form 




Jexp(lnP(A0)L^' 



P{w) = /iexp(-^|w|) 




in which a hyperprior is used given by 
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P(a) = T(a\a 9 b) ^ b a a a - l e" ba /r(a) 
where T(a) is the Gamma function. This distribution has the useful properties 

The marginal distribution of w is then obtained by integrating over a. 
5 The variational framework can be rendered tractable by working not directly with 

the marginal distribution P(w), but instead leaving hierarchical conjugate form explicit, 
and introduce a factorial representation given by Q(w 9 a) = Q(w)Q(a) . This is also 
advantageous because it becomes possible to evaluate the lower bound L as a closed- form 
analytic expression. This is useful for monitoring the convergence of the iterative 

10 optimization, and also for checking the accuracy of an implementation that can be 

embodied in software, such as one or more computer programs running on a computer as 
has been described in conjunction with FIG. 1 in a previous section of the detailed 
description, since updates to the variational distributions cannot lead to a decrease of L. 
It can also be used to compare models, without resorting to a separate validation set, 

15 since it represents an approximation to the model evidence. 

Variational Relevance Vector Regression 

In this section of the detailed description, relevance vector regression, used in 
conjunction with continuous data sets, is described, in accordance with one embodiment 
of the invention, utilizing a variational approach. Thus, the Relevance Vector Machine 
20 (RVM) is augmented by the introduction of hyperpriors given by a separate distribution 
for each hyperparameter a m of the form P(a m ) = T(a m \a,b). Similarly, a prior is 

introduced over the inverse noise variance x given by P{r) = T(r \ c, d). Broad 
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hyperpriors are obtained by setting a = b = c = d =10 in one embodiment. Together 
with the likelihood function 

P(r|X 5 w ? r) = nP(^|x n w 5 r) 

and the weight prior 
5 P{^\a) = f[N(w m \0,a m l ) 

m=0 

a complete probabilistic specification of the model is achieved. 

Next, a factorial approximation is considered to the posterior distribution 
P(yv,a 9 r\X 9 T) given by £>(w,ar,r) = Q^)Q a {a)Q T (r) . Due to the conjugacy 
properties of the chosen distributions, the general solution 

exp(lnP(A£)). 



10 0W) = - 



{exp(lnP(A0)L^ 

can be evaluated analytically, giving 

a(w) = iV(w|// w ,2 w ) 

Q t (t) = t(t\£J) 



m=0 



1 5 where 



S_ = 



diag(^ m ) + (r) 
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a m =a+l/2 b m =b + (w 2 m )/2 
c = c + (N + l)/2 

d-d + \±tl- <w) T + l -±£ (ww T ) A 

b=1 n=l n=l 

5 The requirement moments are evaluated using 

(w) = A w 

(\na m ) = <p(a m )-\nb m 
10 (r) = c/d 

(lnr) = <p(c)-\nd 

where the \|/ function is defined by 

<2?(a) =— lnr(a) 
ry } da W 

The full predictive distribution P(t \ x,X,T) is given by 
15 P(t\x,X,T)= jjp(t\x,yv,r)p(yv,T\X,T)dY/ dr 

In the variational framework, the true posterior P(w, r\X,T) is replaced by the 
variational approximation g w (w) • Q t (t) . Integration over both w and t is intractable. 
However, as the number of data points increases, the distribution of x becomes tightly 
concentrated around its mean value. To see this, it is noted that the variance of x is given 
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by l^ 2 ^ - (t} 2 = c I d 2 ~ 0(1 1 N) for large N. Thus, the predictive distribution can be 
approximated by using 

P(t\x,X,T)= |P^|x ? w ? (r))e w (w)rfw 

which is the convolution of two Gaussian distributions. Using 
5 P(t |x,w ? r) = N(t \y (x,w) , r' 1 ) 

and 

g w (w) = iV(w|// w ,2 w ) 

the following is then obtained 

P(t\x,X,T) = N(t\Sj{x),cr 2 ) 

1 0 where the input-dependent variance is given by 

^(x) = 7 ^ + ^(x) T E^(x) 

W 

The lower bound L can also be evaluated, taking the form 

£ = (lnP(r|X,w,r)) + (lnP(w|tf)) 

+ (lnP(^)) + (lnP(r)>-(lne w (w)) 
-(]nQ a {a))-{\nQ r (r)) 

in which 

15 (lnP(r|X,w4 = ^^ 

^ £ ^ L n=l n=l n=l J 

{i„p(»H)=-^±ih(2.)4i<ta a „)-ii( a .>(^) 

' Z Z m=0 Z m=0 

(lnP(4 = (iV + l)aln& + (a-l)X(lna w )-fe2k)-(^ + l)lnr(a) 

m=0 m=0 
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(lnP(r)) = cln d + (c -l)(lnr) - d (r) - lnr(c) 



(]nQ yf ) = (N + 1)(1 + ln(2*))/2 + ln|l w |/2 



(lna) = tKln^+(a m -l)(b^)-6 m (« m )-br(5 m )} 



m=0 



(ln2 r ) = clnJ+(c-l)(lnr>-J{r)-lnr(c) 



5 Variational Relevance Vector Classification 

In this section of the detailed description, relevance vector regression, used in 
conjunction with discrete data sets, is described, in accordance with one embodiment of 
the invention, utilizing a variational approach. Classification is somewhat more complex 
than regression, because a fully conjugate hierarchical structure is not present. 
1 0 Resolution of this difficulty is now described. The log marginal probability of the target 
data, given the input data, can be written 



As before, a factorized variational posterior of the form £> w (w) • Q a (a) is introduced, and 
the following lower bound on the log marginal probability is obtained 



However, the right-hand side of the above is intractable. Therefore, a further bound is 
introduced using the inequality 



\nP(T\x) = \n jjp(T\X,w)p(w\a)P(a)dw da 



\rvP{T\x)> \\Q^)Q a {a) 



15 
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Z £ / \ 

><x(<?)exp A(s){^z 2 -s 2 J 



where z = (2t- l)y and = (1 /4£) tanh(^7 2). Here, ^ is a variational parameter, 
such that equality is achieved for g = z. Thus, 



P(r|X,w)>F(r,X,w^) = ^ 

«=i V ^ 



5 where z n = (2t n - l)wV„- Substituting into 

]nP(T\x)> \jQ w (yv)Q a (a) 

1 <2 w (w)Q» J 

and noting that P(T \ X,w)/ F(T,X,w,4) > 1 implies \nP(T \ X,w)/ F(T,X,w,g) > 0 , 

a lower bound on the original lower bound is obtained, 

\nP(T\x) > L = jjdw daQ v (w)g a («) 
f F(r,X,w)P(w|g)P(#) ' 

1 0 The right-hand side of the above equation is now optimized with respect to 

Q w (w) and Q a (a) , as well as with respect to the parameters # = The variational 
optimization for Q w (w) yields a normal distribution of the form 

G w (w) = tf(w|m,S) 



S = 



A + 2|>(* B )^ T 



15 



m= i s fZ( 2 ^- 1 )^ 
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where A = diag{# m ) . Similarly, variational optimization of Q a (a) yields a product of 



Gamma distributions of the form 



10 



15 



m=0 



a = a + — 
2 



Finally, maximizing 

lnP(T\x) >L= \\dw daQ w {w)Q a (a) 
' F(T,X,w)P(w\a)P(a) 



ln< 



with respect to the variational parameters % n gives re-estimation equations of the form 

The lower bound given by the right-hand side of 

\nP(T\x)>L= l\dw daQ„(w)Q a (a) 

ln< 



F(T,X,w)P(w\a)P(a) 



Q m (*)Q.(a) 

can also be evaluated, as 

L = (lnF) + (lnP(w|ar)) + (hP(ff)) - (infi. (w)) - (Q a (a)) 



where 



( m F> = £{ln^>±(2^-^ 

(toP(w|4 = 4i(^)(^)4t(ln«J-^ln(2^) 
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(in P(a)) = £ [-ba/b + (a - 1) (q> (a) - In b ) + a In b - In r (a)} 

m=0 

-(ln&(w)) = ^(l + ln2*) + ±ln|S| 

-(lna^^^J-f^-O^aJ-ln^+^ + lnr^)} 

Predictions from the trained model for new inputs can be obtained by substituting 
5 the posterior mean weights into 

P(t\x,v) = cr(yy[l-a(y)y 

to give the predictive distribution in the form 

P(f|*,(w)). 

Furthermore, a more accurate estimate can take into account the weight uncertainty by 
10 marginalizing over the posterior distribution of the weights. Using the variational result 
g w (w) for the posterior distribution leads to convolution of a sigmoid with a Gaussian, 

which is intractable. From symmetry, however, such a marginalization does not change 
the location of the p - 0.5 decision surface. 

Methods 

15 In this section of the detailed description, methods for implementing a variational 

relevance vector machine, according to varying embodiments of the invention, are 
presented. The variational relevance vector machine allows a user to generate a simpler 
(less complex) model as compared to other techniques, such as a support vector machine, 
using a variational approach. The computer-implemented methods are desirably realized 

20 at least in part as one or more programs running on a computer - that is, as a program 



25 



executed from a computer-readable medium such as a memory by a processor of a 
computer. The programs are desirably storable on a machine-readable medium such as a 
floppy disk or a CD-ROM, for distribution and installation and execution on another 
computer. The program or programs can be a part of a computer system or a computer, 
5 such as that described in conjunction with FIG. 1 in a previous section of the detailed 
description. It is noted that the description of the methods in this section implicitly rely 
on the description of the various variational relevance vector machine techniques as has 
been described in the preceding sections of the detailed description, as well as other of 
the preceding sections of the detailed description. 
Q 10 Referring to FIG. 4, a flowchart of a method 399 according to an embodiment of 

^ ™i 

!; ["! the invention is shown. In 400, a data set to be modeled is input. The data set desirably 

ry 

5~ represents a real-world event or events - for example, the stories read by users of a web 

rrj site, the television shows watched by viewers, etc. Inputting the data set in one 

O embodiment also includes selecting a model for the data set, such as a generalized linear 

Q 1 5 model as has been described - that is, a form for the ultimate model for the data set can 

□ be selected. The invention is not limited to the manner by which input is conducted in 

.; I - 
z. -.^ 

400. For example, the data set may be input by being received from a communications 
device, from a storage device, etc. 

In 402, an initial set of hyperparameters is selected. Ultimately, the set of 
20 hyperparameters is used for determining a prior distribution for the data set, for modeling 
thereof, using a variational approach. The prior distribution is approximated by a product 
of the distribution of the set of hyperparameters, the distribution of a set of weights, and 
the distribution of a set of predetermined additional parameters, which in the case of a 
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continuous data set for regression are parameters accounting for noise, and which in the 
case of a discrete data set for classification are parameters accounting for a (second) 
lower bound. 

Thus, 404, 406, 408 or 41 0, and 412 represent the iterative cycle of the variational 

5 approach for a relevance vector machine. In 404, the distribution for the set of weights is 
updated, and then in 406, the distribution of the set of hyperparameters (the a 
parameters) is updated. Next, depending on whether regression or classification is being 
performed, either the distribution for the parameters taking into account noise (the % 
parameters), in the former case, or the distribution for the parameters taking into account 

10 the (second) lower bound (the £ parameters). If the predetermined convergence criterion 
(as alluded to in the preceding sections of the detailed description) has not been achieved, 
then the method iterates again, proceeding from 412 back to 404. Otherwise, the method 
proceeds from 412 to 414. 

In 414, at least the posterior distribution, as has been determined, is output. The 

15 complete model, including the distribution, may be output as well in one embodiment. 
The invention is not limited to the manner by which output is accomplished. For 
example, in one embodiment, output can be to a further analysis program or software 
component, that allows for analysis and conclusions to be drawn. As another example, 
output can be displayed on a displayed device, or printed to a printer, etc. As a third 

20 example, output can be stored on a storage device, for later further analysis program or 
software component. 
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Conclusion 

Although specific embodiments have been illustrated and described herein, it will 
be appreciated by those of ordinary skill in the art that any arrangement which is 
calculated to achieve the same purpose may be substituted for the specific embodiments 
shown. This application is intended to cover any adaptations or variations of the present 
invention. Therefore, it is manifestly intended that this invention be limited only by the 
following claims and equivalents thereof. 
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We claim: 



1 . A computer-implemented method for modeling a data set comprising: 

selecting an initial set of hyperparameters for determining a prior distribution for the 

data set for modeling thereof, the prior distribution approximated by a product of a 
5 distribution of the set of hyperparameters, a distribution of a set of weights, and a 

distribution of a set of predetermined additional parameters; and, 

iteratively updating the distribution of the set of weights, the distribution of the set of 

hyperparameters, and the distribution of the set of predetermined additional parameters 

until a predetermined convergence criterion has been reached, 
1 0 such that the product of the distribution of the set of hyperparameters, the distribution 

of the set of weights, and the distribution of the set of predetermined additional 

parameters as have been iteratively updated approximates the posterior distribution for 

modeling of the data set. 

2. The method of claim 1 , further comprising initially inputting the data set to be 
15 modeled. 

3 . The method of claim 1 , further comprising outputting the posterior distribution as 
approximated by the product of the distribution of the set of hyperparameters, the set of 
weights, and the set of predetermined additional parameters. 
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4. The method of claim 1, wherein the data set comprises a continuous data set, such 
that the set of predetermined additional parameters comprises a set of parameters 
accounting for noise. 

5. The method of claim 1, wherein the data set comprises a discrete data set, such that 
5 the set of predetermined additional parameters comprises a set of parameters accounting 

for a lower bound. 

6. A computer-implemented method comprising: 
inputting a data set to be modeled; 

determining a relevance vector learning machine via a variational approach to obtain 
1 0 a posterior distribution for the data set; and, 

outputting at least the posterior distribution for the data set. 

7. The method of claim 6, wherein determining the relevance vector learning machine 
comprises: 

selecting an initial set of hyperparameters for determining the prior distribution for 
1 5 the data set, the prior distribution approximated by a product of a distribution of the set of 
hyperparameters, a distribution of a set of weights, and a distribution of a set of 
predetermined additional parameters; and, 

iteratively updating the distribution of the set of weights, the distribution of the set of 
hyperparameters, and the distribution of the set of predetermined additional parameters 
20 until a predetermined convergence criterion has been reached. 
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8. The method of claim 6, wherein the data set comprises a discrete data set, such that 
the set of predetermined additional parameters comprises a set of parameters accounting 
for a lower bound. 

9. The method of claim 6, wherein the data set comprises a continuous data set, such 
5 that the set of predetermined additional parameters comprises a set of parameters 

accounting for noise. 

10. A machine-readable medium having instructions stored thereon for execution by a 
processor to perform a method for modeling a continuous data set comprising: 

selecting an initial set of hyperparameters for determining a prior distribution for the 
10 continuous data set for modeling thereof, the prior distribution approximated by a product 
of a distribution of the set of hyperparameters, a distribution of a set of weights, and a 
distribution of a set of parameters accounting for noise; and, 

iteratively updating the distribution of the set of weights, the distribution of the set of 
hyperparameters, and the distribution of the set of parameters accounting for noise until a 
1 5 predetermined convergence criterion has been reached, 

such that the product of the distribution of the set of hyperparameters, the distribution 
of the set of weights, and the distribution of the set of parameters accounting for noise as 
have been iteratively updated approximates the posterior distribution for modeling of the 
continuous data set. 

20 11. The medium of claim 10, wherein the method further comprises initially inputting the 
continuous data set to be modeled. 
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12. The medium of claim 10, wherein the method further comprises further comprising 
outputting the posterior distribution as approximated by the product of the distribution of 
the set of hyperparameters, the set of weights, and the set of parameters accounting for 
noise. 

5 1 3 . A machine-readable medium having instructions stored thereon for execution by a 
processor to perform a method for modeling a discrete data set comprising: 

selecting an initial set of hyperparameters for determining a prior distribution for the 
discrete data set for modeling thereof, the prior distribution approximated by a product of 
a distribution of the set of hyperparameters, a distribution of a set of weights, and a 
10 distribution of a set of parameters accounting for a lower bound; and, 

iteratively updating the distribution of the set of weights, the distribution of the set of 
hyperparameters, and the distribution of the set of parameters accounting for a lower 
bound until a predetermined convergence criterion has been reached, 

such that the product of the distribution of the set of hyperparameters, the distribution 
15 of the set of weights, and the distribution of the set of parameters accounting for a lower 
bound as have been iteratively updated approximates the posterior distribution for 
modeling of the discrete data set. 

14. The medium of claim 10, wherein the method further comprises initially inputting the 
discrete data set to be modeled. 

20 15. The medium of claim 10, wherein the method further comprises further comprising 
outputting the posterior distribution as approximated by the product of the distribution of 
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the set of hyperparameters, the set of weights, and the set of parameters accounting for a 
lower bound. 



16. A machine-readable medium having instructions stored thereon for execution by a 
processor to perform a method comprising: 

5 inputting a data set to be modeled; 

determining a relevance vector learning machine via a variational approach to obtain 
a posterior distribution for the data set; and, 

outputting at least the posterior distribution for the data set. 

17. The medium of claim 16, wherein determining the relevance vector learning machine 
10 comprises: 

selecting an initial set of hyperparameters for determining the prior distribution for 
the data set, the prior distribution approximated by a product of a distribution of the set of 
hyperparameters, a distribution of a set of weights, and a distribution of a set of 
predetermined additional parameters; and, 
15 iteratively updating the distribution of the set of weights, the distribution of the set of 
hyperparameters, and the distribution of the set of predetermined additional parameters 
until a predetermined convergence criterion has been reached. 

18. The medium of claim 16, wherein the data set comprises a discrete data set, such that 
the set of predetermined additional parameters comprises a set of parameters accounting 

20 for a lower bound. 
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19. The medium of claim 16, wherein the data set comprises a continuous data set, such 
that the set of predetermined additional parameters comprises a set of parameters 
accounting for noise. 



34 



ABSTRACT OF THE DISCLOSURE 

A variational Relevance Vector Machine (RVM) is disclosed. The RVM is a 
probabilistic basis model. Sparsity is achieved through a Bayesian treatment, where a 
prior is introduced over the weights governed by a set of what are referred to as 

5 hyperparameters - one such hyperparameter associated with each weight. An 

approximation to the joint posterior distribution over weights and hyperparameters is 
iteratively estimated from the data. The posterior distribution of many of the weights is 
sharply peaked around zero, in practice. The variational RVM utilizes a variational 
approach to solve the model, in particular using product approximations to obtain the 

10 posterior distribution. 

I hereby certify that this is being deposited with the United States 
Postal Service "Express Mail Post Office to addressee" service 
under 37 CFR § 1 .10 in an envelope addressed to The Assistant 
Commissioner for Patents, Washington, DC 20231, on Mar 7, 



2000 , by Sherry R. Drvia , and having "express mail" mailing label 
no. EK530297161US . 




35 




LU ^ 
c£ ad 
< O 

*S 

O LU 

O 2 



o 

CM 



oo 



oo 



lu \r 

> Q 
< 



V LU 

8k 



to 



__j o 

£££ 

CO | — 



CM. 



O 

CO 
CO 
LU 

o 
o 
a: 

Q_ 



CO 
=) 

LU 
h— 

CO 

CO 



V 



LU 

— I o 

< LU < 
O > LL 

F cE °^ 



o 



°°k\ 



CM 



o 



LU 

, O 

i — LU <f 

lu ^ > S: 

2: ^ £ 

O Q 9k LU 

< H- 



00 



1' 



CM 



M^ LU 

CO o 

A ty < 

Q > LL_ 

Q 5 



CM_^_ 



□ 
□ 



CD 



O! 



LU 

I— 

CO 
CO 1 



o 

Of 



CO 

o 

OQ 



LO 
CO. 



CO 
CO 



0 








1 


LU 




h- 




CO 


LU 


> 


CL 


CO 


O 





1. 



OO 
CO. 



LOi 
CM I 



II 
<°- 



S CO 
CC < LLj 
1X1 CC =^ 

— o 2 

. 06 
o 2 Q 



8S 

CL 



or 
o 



LU 



XL 

< 

LU 
Q 



3 







LU LU 




1— I — 




O -O 








LU ^> 




cr: 0 




0 





00 



CNi 



a> 
cm^ 



oo r 



ION 


MS 


r— 




< 




O 






0 


APP 


PR 



\o 

^3- 



o 
o 

Q_ 



< 
Q 



< 

CD 
O 

a: 
□_ 



CO 
LU 



o 



CO 



< 
0 



CD 



h- LU 

CO 
LU > 
Q_ CO 
O 



FIG 2 



INPUT 



202, 




OUTPUT 
204 



FIG 3 



LEARNING 
MACHINE 



300, 



p(c|x) 
304, 



302. 



399- 




400 



FIG 4 



CHOOSE 
HYPERPARAMETERS 



402 



UPDATE WEIGHT 
DISTRIBUTION 



404 



UPDATE ALPHA 
DISTRIBUTION 



406 



REGRESSION 



CLASSIFICATION 



UPDATE TAU 
DISTRIBUTION 



UPDATE XI 
PARAMETERS 



Lr\ 410 




PATENT APPLICATION 



Declaration and power of attorney 
for patent application 



ATTORNEY DOCKET NO. 1018.093US1 



As a below named inventor, I hereby declare that: 

My residence/post office address and citizenship are as stated below next to my name; 

I believe I am the original, first and sole inventor (if only one name is listed below) or an original, first and joint inventor (if plural names 
are listed below) of the subject matter which is claimed and for which a patent is sought on the invention entitled: 
Variational Relevance Vector Machine 

the specification of which is attached hereto unless the following box is checked: 

( ) was filed on as US Application Serial No. or PCT International Application 

Number and was amended on (if applicable). 

I hereby state that I have reviewed and understood the contents of the above-identified specification, including the claims, as amended by 
any amendments) referred to above. 1 acknowledge the duty to disclose all information which is material to patentability as defined in 
37CFRL56. 

Foreign Application^) and/or Claim of Foreign Priority 

I hereby claim foreign priority benefits under Title 35, United States Code Section 119 of any foreign application(s) for patent or mventor(s) certificate listed below and 



COUNTRY 


APPLICATION NUMBER 


DATE FILED 


PRIORITY CLAIMED UNDER 35 U.S.C. 119 








YES: NO- 








YES: NO: 



Provisional Application 

I hereby claim the benefit under Title 35, United States Code Section 1 19(e) of any United States provisional application(s) listed below: 



APPLICATION SERIAL NUMBER 


FILING DATE 











; $.S. Priority Claim 

^hereby claim the benefit under Title 35, United States Code, Section 120 of any United States application(s) listed below and, insofar as the subject matter of each of the 
; ieriaims of this application is not disclosed m the prior United States application m the manner provided by the first paragraph of Title 35, United States Code Section 112, 
['^acknowledge the duty to disclose material information as defined in Title 37, Code of Federal Regulations, Section 1.56(a) which occurred between the filing date of the 
i!§rior application and the national or PCT international filing date of this application: 



j APPLICATION SERIAL NUMBER 


FILING DATE 


ST ATUS (paten te d/pend mg/aband on e d) 





















[ JOWER OF ATTORNEY: 

I'^ls a named inventor, I hereby appoint the following attorney(s) and/or agent(s) listed below to prosecute this application and transact all business in the Patent and 
trademark Office connected therewith. 



Michael A. Dryja, Reg. No. 39662 



Katie E. Sako, Reg. No. 32628 



Daniel D. Crouse, Reg. No. 32022 



Send Correspondence tou , 


Direct Telephone Calls To: 


Michael A. Dryja 


Michael A. Dryja 


Law Offices of Michael Dryja 


425-427-5094 


704 228th Avenue NE 




PMB694 




Redmond, WA 98053 





I hereby declare that all statements made herein of my own knowledge are true and that all statements made on information and belief are believed to be true, and further 
that these statements were made with the knowledge that willful false statements and the like so made are punishable by fine or imprisonment, or both, under Section 
1001 of Title 18 of the United States Code and that such willful false statements may jeopardize the validity of the application or any patent issued thereon. 



Full Name of Inventor: Michael Tipping 

Residence: 1 The Brambles, Girton, Cambridge, England CB3 ONY United Kingdom 
Post Office Address: Same 



Citizenship: British 



Inventor's Signature 



Date 



Page 1 of 2 



DECLARATION AND POWER OF ATTORNEY 



MR PATENT APPUCATIO^ntiBHed) ATTORNEY DOCKET NO. 1018.093US1 



Full Name of Inventor: Christop her Bishnp 

Citizenship: British 

Residence: ^LYndewodeKoark Ca mbrid^ En p i an d CB1 2HN VriteA Ki nj ^ 
Post Office Address: Same 



Inventor's Signature 



Date 



I hereby certify that this is being deposited with the United States 
Postal Service "Express Mail Post Office to addressee" service 
under 37 CFR § 1.10 in an envelope addressed to The Assistant 
Commissioner for Patents, Washington, DC 20231, on Mar 7, 
2000 , by Sherry R, Drvia , and having "express maii" mailing label 
no. EK530297161US . 




Signature of Starry R. Dryja 



Page 2 of 2 



